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The star ( Ophiuchi is one of the brightest massive stars in the northern hemisphere and was intensively studied in various 
wavelength domains. The currently available observational material suggests that certain observed phenomena are related 
to the presence of a magnetic field. We acquired spectropolarimetric observations of £ Oph with FORS 1 mounted on 
the 8-m Kueyen telescope of the VLT to investigate if a magnetic field is indeed present in this star. Using all available 
absorption lines, we detect a mean longitudinal magnetic field {B Z ) M = 141 ± 45 G, confirming the magnetic nature 
of this star. We review the X-ray properties of £ Oph with the aim to understand whether the X-ray emission of £ Oph is 
dominated by magnetic or by wind instability processes. 
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1 Introduction 

During the last years a gradually increasing number of O, 
early B-type, and WR stars have been investigated for mag- 
netic fields, and as a result, about a dozen magnetic O-type 
stars are presently known (e.g., Hubrig et al. 2008; Martins 
et al. 120101 Hubrig et al. 1201 la) . The recent detections of 
magnetic fields in massive stars generate a strong motiva- 
tion to study the correlations between evolutionary state, ro- 
tation velocity, and surface composition, and to understand 
the origin and the role of magnetic fields in massive stars. 

The star C Ophiuchi (=HD 149757) of spectral type 
09.5V is a well-known rapidly rotating runaway star with 
extremely interesting characteristics. It undergoes episodic 
mass loss seen as emission in Ha, and it is possible that 
, it rotates with almost break-up velocity with vsini = 
400kms _1 (Kambe et al. 11993) . Various studies indicate 
different types of spectral and photometric variability. The 
UV resonance lines show multiple discrete absorption com- 
ponents (DAC) in the UV (e.g. Howarth et al. 11984) and 
strong line profile variations in optical spectra reconciled 
with traveling sectorial modes of high degree (e.g. Reid et 
al. 11993) . Highly precise MOST (Microvariability and Os- 
cillations of Stars) satellite photometry in 2004 has yielded 
at least a dozen significant oscillation frequencies between 
1 and 10 cycles/day, hinting at a behaviour similar to 
j3 Cephei-type stars (Walker et al. 2005). No unambiguous 
rotation period could be identified in spectroscopic and pho- 

* Based on observations obtained at the European Southern Observa- 
tory, Paranal, Chile (ESO programme 081.C-0410(A)). 
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tometric observations, although Balona & Kambe (11999) fa- 
vored a period in the region of 1 cycle/day. 

COph is also well-known for its variability in the X- 
ray band. Oskinova et al. (12001) studied the ASCA obser- 
vations of COph that covered just more than the expected 
rotation period of the star. A clearly detected periodic X- 
ray flux variability with ^20% amplitude in the ASCA pass- 
band (0.5-10keV) was reported. A period of 0.77 d was de- 
tected and a possible connection with the recurrence time 
(0.875 d±0. 167 d) of the DACs in UV spectra of the star was 
discussed. The DACs in the spectra of O stars are commonly 
explained by large-scale structures in the stellar wind, mod- 
ulated by rotation and possibly related to a surface mag- 
netic field (Cranmer & Owocki 1996). Waldron et al. (in 
preparation, private communication) found that SUZAKU 
data on ( Oph suggest a period of ^0.98 d that is consis- 
tent but slightly larger than the X-ray periodicities found in 
ASCA data (Oskinova et al. [200TT ) and in Chandra HETGS 
data (Waldron |200"5l In addition, the HETGS data appear to 
indicate an additional cyclic period of ^0.33 d in the hard 
X-ray band (>1. 2 keV). 

The results of our previous studies seem to indicate that 
the presence of a magnetic field is more frequently detected 
in candidate runaway stars than in stars belonging to clusters 
or associations (Hubrig et al. 1201 lbi Hubrig et al. 1201 la) . 
The currently best available astrometric, spectroscopic, and 
photometric data were used to calculate the kinematical sta- 
tus of magnetic O-type stars with previously unknown space 
velocities. The results suggest that most of the magnetic O- 
type stars can be considered as candidate runaway stars. 
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The available observational material suggests that ( Oph 
is a main sequence single star in the field with runaway char- 
acteristics. Usually, to explain the origin of massive stars in 
the field, two mechanisms are discussed in the literature. In 
one scenario, close multibody interactions in a dense cluster 
cause one or more stars to be scattered out of the region (e.g. 
Leonard & Duncan 1990). For this mechanism, runaways 
are ejected in dynamical three- or four-body interactions. 
An alternative mechanism involves a supernova (SN) explo- 
sion within a close binary, ejecting the secondary due to the 
conservation of momentum (Zwicky |l9571 Blaauw [l9611 ). 
Blaauw (1952) suggested the origin of COph in the Scor- 
pius OB2 association due to its proper motion vector, which 
points away from the association. More recently, Hoogerw- 
erf et al. (2001) suggested that the star gained it space veloc- 
ity of ~30kms _1 in a supernova explosion within a close 
binary in Upper Scorpius about 1-2 Myr ago. The authors 
identified PSR B 1929+ 10 as an associated pulsar with a 
characteristic age of ^3 Myr, consistent with the kinematic 
age of COph within the uncertainties. Tetzlaff et al. (2010) 
reinvestigated the scenario of a binary SN in Upper Scorpius 
involving COph and PSR B 1929+ 10 and concluded that it 
is very likely that both objects were ejected during the same 
supernova event. In their work, the considered association 
age range implies that the progenitor star of the produced 
neutron star had a spectral type between 06/07 and 09 
with a mass range from 18 to 37 M®. The X-ray emission 
of the pulsar seems to be dominated by non-thermal radia- 
tion processes (e.g. Becker et al. 2006). An arc-like nebula 
surrounding PSR B 1929+ 10 and extending up to 10" was 
identified in Chandra data and interpreted as a bow-shock 
nebula (Hui & Becker 2008 ). The estimated magnetic field 
strength in the shocked region accounts for ^75 pG, while 
the typical magnetic field strength in the ISM is about 2- 
6 /iG. 

The presence of a bow-shock nebula has also been de- 
tected for C Oph. FigureQ]shows an image based on archival 
Spitzer IRAC maps (AOR 17774848). Recently, Kobulnicky 
et al. ( 120101 ) analyzed a sample of bow shocks around mas- 
sive stars in Cygnus-X. They used the analytical descrip- 
tion of momentum-driven bow shocks and dust/polycyclic 
aromatic hydrocarbon emission models to estimate stellar 
mass loss rates from the observed properties of the bow 
shocks. It was found that mass-loss rates in the range be- 
tween 10~ 7 M Q yr -1 and a few times 10 _6 M Q yr^ 1 are 
required to generate the bow shocks around typical B2V - 
05V type stars. 

The mass-loss rate M of this star was empirically ob- 
tained from different diagnostics by a number of authors. 
Repolust et al. (2004 ) fitted the Ha photospheric absorp- 
tion line and derived the upper limit on the C Oph mass-loss 
rate as 1.8 x 10~ 7 A/ o yr -1 . Fullerton et al. (120061) deter- 
mined the radio-based mass-loss rate of COph as 1.1 x 
10 _7 A/ Q yr _1 . The mass-loss rates determined from ra- 
dio depend on the square of the density since the physical 
mechanism responsible for the radio emission is free-free 




Fig. 1 (online colour at: www.an-journal.org) Combined 
IR Spitzer IRAC (3.6 /im blue, 4.5 pm green, 8.0 /im red) 
image of the bow shock around the runaway star COph. 
Archival data have been used. Galactic coordinates are 
shown. The image size is ~ 36' x 31'. 



emission. On the other hand, Fullerton et al. derive a much 
smaller mass-loss rate from fitting the UV P V resonance 
doublet, the product of the mass-loss rate and the ion frac- 
tion of P+ 4 being only Mg(P+ 4 )<1.3 x lO^A/syr" 1 
with q(P +4 )<l. The mass-loss rates derived from fitting 
the wind profiles of UV resonance lines depend linearly on 
the density. To resolve this discordance in mass-loss deter- 
minations based on p 2 - and p-diagnostics, Fullerton et al. 
suggest that the winds are strongly clumped with a volume 
filling factor of — 10 3 -10 5 . Marcolino et al. (2009b an- 
alyzed optical and UV spectra of C Oph among their sam- 
ple of O-type dwarfs. They derive an upper limit on the 
mass-loss rate of C Oph as 1.6 x 10 _9 Mq yr _1 if the wind 
was smooth. This value agrees with the Mq(P +i ) value ob- 
tained by Fullerton et al. (120061 ) . 

Using the example of the O-type supergiant CPuppis, 
Oskinova et al. (2007) demonstrated that the discordance of 
mass-loss rates found by Fullerton et al. can be overcome 
by accounting for stellar wind porosity (see also Sundqvist 
et al. 120 101 ). It was found for the 051a star C Puppis that only 
a moderate reduction of the mass-loss rate by a factor of 2-3 
(compared to the smooth wind models) is required to repro- 
duce both Ha and P V lines. If this result holds also for non- 
supergiant O type stars, the mass-loss rate of C Oph is only 
a few times lower compared to the radio-based mass-loss 
determined by Fullerton et al., i.e. ~10 -7 Mq yr^ 1 . Impor- 
tantly, this mass-loss rate is in agreement with values that 
are required to produce bow shocks around O stars (Kobul- 
nicky et al. 1201 01 ). 

An additional aspect, which may hint at the presence 
of a magnetic field in runaway stars, is that a number of 
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individual abundance studies indicate nitrogen enrichment 
in the atmospheres of runaway stars (e.g. Boyajian et al. 
2005). Nitrogen enrichment was found in COph by Villa- 
mariz & Herrero (2005). Recent NLTE abundance analy- 
ses (e.g., Morel et al. 120081 Hunter et al. 120081 1 suggest that 
slow rotators have peculiar chemical enrichment such as ni- 
trogen excess or boron depletion, and these peculiarities are 
linked to the presence of a magnetic field. On the other hand, 
Hubrig et al. (1201 lei showed that some magnetic massive 
stars previously assumed to be slow rotators, are in fact fast 
rotators, but are viewed close to their rotation poles. 

To test the magnetic nature of this particularly interest- 
ing rapidly rotating runaway star, we acquired spectropo- 
larimetric observations with the low-resolution spectropo- 
larimeter FORS 1 at the VLT. In this work we report the 
first detection of a magnetic field in this star. 

2 Magnetic field measurements 

Spectropolarimetric observations with FORS 1 (Appen- 
zeller et al. [19981 were obtained on 2008 May 23 
(MJD 54609.34) using grism 600B and a slit width of 0'/4 
to achieve a spectral resolving power of R w 2000. The 
use of the mosaic detector made of blue optimized E2V 
chips and a pixel size of 15 /im allowed us to cover a large 
spectral range, from 3250 to 6215 A, which includes all 
hydrogen Balmer lines from H/3 to the Balmer jump. Six 
continuous series of two exposures with durations between 
0.3 and 3 sec were taken at two retarder waveplate setups 
(a = +45° and —45°). For all but the first exposure pairs 
we achieved a signal-to-noise ratio between 1000 and 1500. 
More details on the observing technique with FORS 1 can 
be found elsewhere (e.g., Hubrig et al. I2004al I2004bl and 
references therein). The mean longitudinal magnetic field, 
(B z ), was derived using 

V 
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where V is the Stokes parameter that measures the circu- 
lar polarisation, / is the intensity in the unpolarised spec- 
trum, g e g is the effective Lande factor, e is the electron 
charge, A is the wavelength, m e the electron mass, c the 
speed of light, d//dA is the derivative of Stokes /, and (B z ) 
is the mean longitudinal magnetic field. The longitudinal 
magnetic field was measured in two ways: using only the 
absorption hydrogen Balmer lines or using the entire spec- 
trum including all available absorption lines. We obtain for 
the mean longitudinal magnetic field using all available ab- 
sorption lines (B z ) all — 141 ± 45 G and for the mean lon- 
gitudinal magnetic field using the hydrogen Balmer lines 
(-Bz)hyd = 123 ± 54G. Our detection using the entire spec- 
trum has a significance of 3.1a, determined from the for- 
mal uncertainties we derive. In the Stokes V spectra dis- 
tinct Zeeman signatures are well visible at the position of 
hydrogen and metal lines. In Fig.|2]we display Stokes / and 
V spectra in the spectral regions around H/3 and the Na I 
doublet. 
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Fig. 2 Stokes / and V spectra of £ Oph in the spectral 
regions around H/3 (top) and the Na I doublet (bottom). 

In Fig. [3] we present time series of Stokes / spectra cor- 
responding to our data set of six sub-exposures in the region 
around He II 4686 A and He I 4713 A. Although the time 
lapse between the observations of the first pair and the last 
pair is only 13 minutes, some small line profile variations in 
the He I line are already detectable at such short time scales. 

The only other measurements of the magnetic field in 
COph have been presented by Schnerr et al. (2008) who 
used the MuSiCoS spectropolarimeter to derive Stokes / 
and V spectra with the Least Square Deconvolution method. 
No longitudinal magnetic field was detected in this star at 
a level more than 3a. However, the measurement errors in 
their observations were in the range from 700 G to more 
than 3 kG. 

Massive stars usually end their evolution with a final su- 
pernova explosion, producing neutron stars or black holes. 
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Fig. 3 Six Stokes / spectra corresponding to six sub- 
exposures in the region around He II 4686 A and He I 
4713 A. Small line profile variations in the He I line are al- 
ready detectable on the time scale of only 13 min. 

The initial masses of these stars range from ~8— 10 M@ to 
100M Q or more, which correspond to spectral types ear- 
lier than about B2. Contrary to the case of Sun-like stars, 
the magnetic fields of stars on the upper main sequence 
(Ap/Bp stars), white dwarfs, and neutron stars are domi- 
nated by large spatial scales and do not change on yearly 
time scales. In each of these classes there is a wide distri- 
bution of magnetic field strengths, but the distribution of 
magnetic fluxes appears to be similar in each class, with 
maxima of <£> max = ttR 2 B ~ 10 27_28 G cm 2 (Reiseneg- 
ger l200lt Ferrario & Wickramasinghe 2005 ), arguing for a 
fossil field whose flux is conserved along the path of stel- 
lar evolution. According to Reisenegger (2009) the mag- 
netic fluxes have possibly been generated on or even before 
the main-sequence stage and then inherited by the compact 
remnants. 

The magnetic field strength of the pulsar PSR B 1929+10 
is 0.5129x 10 12 Gauss. Assuming simple conservation of 
magnetic flux we obtain a field strength of just a few Gauss 
for the more massive pulsar progenitors. For COph, which 
is supposed to be formed in a binary SN, sharing the same 
environment with the SN progenitor, the expected field 
strength would be of the order of 10 G. This value is no- 
tably lower than our current measurement, possibly indicat- 
ing that either the magnetic field of this middle-aged pulsar 
has significantly decayed during the few Myrs after the SN 
explosion or other mechanisms play a role in the generation 
of magnetic fields in O-type stars. 

3 Discussion 

^Oph has been extremely well studied in all wavelength 
ranges, from the X-ray by all major X-ray satellites (with 
the exception of XMM-Newton) to the infrared region with 
Spitzer. In view of the detection of a magnetic field on ( Oph 
reported in this work, we review its X-ray properties with 



the aim to understand whether the X-ray emission of £ Oph 
is dominated by magnetic or wind instability processes. 

Babel & Montmerle ( 119971 ) studied the case of a rotat- 
ing star with a dipole magnetic field sufficiently strong to 
confine stellar wind. The magnetic field locally dominates 
the bulk motion of stellar wind, when the ratio of mag- 
netic to kinetic energy density, B 2 / [lopv 2 > 1, where v 
is the supersonic flow speed. A collision between the wind 
streams from the two hemispheres in the closed magneto- 
sphere leads to a strong shock and X-ray emission. 

MHD simulations in the framework of this magnet- 
ically confined wind shock (MCWS) model were per- 
formed by ud-Doula & Owocki (2002) and Gagne et al. 
(2005). Using their notation, the wind is confined when 
7/„ = (i? 2 i? 2 )(Afu 0O ) _1 > 1. New observations are re- 
quired to establish whether the magnetic field of £ Oph is 
a dipole. However, for the purpose of this discussion, let 
us assume that the field has an average strength of 150 G. 
Using the stellar parameters of £ Oph as inferred by Mar- 
colino etal. (120091 ), we estimate 77* (£ Oph) ~ 10 3 , i.e. the 
magnetic field should dominate the wind motion up to the 
Alfven radius that is located at <10i?». In this case, the X- 
ray emission should mainly originate from the MCWS. 

The MCWS model predicts that the X-ray emitting 
plasma should be located at a few i?» from the photosphere; 
that the X-ray emission lines should be narrow; that the X- 
ray luminosity should be higher and the spectrum harder 
than in non-magnetic stars; that in oblique magnetic rota- 
tors the X-ray emission should be modulated periodically 
as a consequence of the occultation of the hot plasma by a 
cool torus of matter, or by the opaque stellar core. 

The lines of He-like ions observed in X-ray spectra are 
useful to derive the location of the line formation region in 
hot stars because forbidden line emission is depressed by 
ultraviolet pumping. The latter depends on the distance to 
the stellar photosphere (Gabriel & Jordan 1969; Blumen- 
thal et al. 119721 ). The Si XIII line observed in the Chan- 
dra HETGS/MEG spectrum is shown in Fig. [4] prominent 
forbidden line can easily be distinguished in this fugure, 
while normally forbidden lines are strongly suppressed in 
the spectra of OB-type stars. The presence of the forbid- 
den line implies that the line formation region is located far 
from the photosphere, so that the radiative excitation does 
not lead to the depopulation of the corresponding metastable 
energy levels. Waldron & Cassinelli (2008) found that the 
Si XIII line is formed at 1.8 ± 0.7 i?* in ( Oph and that other 
He-like lines are formed even further out in the wind. In- 
terestingly, the strong forbidden Si XIII line is also observed 
in the Chandra spectrum of the magnetic star r Sco. Cohen 
et al. (120051 ) derive a Si XIII line formaiton radius for r Sco 
in the range between 1.1 i?* and 1.5 R*. These radii of line 
formation are smaller than those found in the prototypical 
MCWS model object (^OriC, 3.4 ± 0.8 (Waldron & 
Cassinelli[2008]). 

Oskinova et al. (2006 ) studied the Chandra spectrum of 
£ Oph among other O-type stars. They found that the X-ray 
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Fig. 4 The Si XIII line observed in the spectrum of £ Oph 
(co-added MEG ±1). Vertical dashed lines indicate the rest- 
frame wavelength: Ar - resonant line, Ai - sum of in- 
tercombination lines, Af - forbidden line. The rest-frame 
wavelengths are corrected for the radial velocity taken from 
Hoogerwerf et al. (12001b . 



emission lines in this star are narrow and that the signatures 
of wind absorption on line profiles are weak. Figure|5]shows 
the Fe XVII and Ne X lines as measured by Chandra plotted 
over units of the wind terminal velocity, w oo =1550kms _1 . 
The lines are only slightly broadened, if at all. 

The X-ray luminosity of C Oph, Lx = 1.2 X 10 31 erg s -1 
and the ratio Lx/^boi = 4 x 10~ 8 are quite usual among 
late type OV stars (Oskinova et al. [2006). Adopting the 
mass-loss rate from COph as <1.8 x lO _7 M0yr _1 , Os- 
kinova et al. (2006) noticed that in COph the ratio of 
X-ray to wind mechanical luminosity L mcc h (Mw 00 2 /2), 
^x/^mech^8.5 x 10 -5 , is a few times higher than in other 
single O-type stars. This may be related to the lower wind 
opacity in C Oph, or it may hint at some additional mecha- 
nism of X-ray generation besides the intrinsic wind shocks. 

From their analysis of Chandra spectra, Zhekov& Palla 
(2007) derived the differential emission measure (DEM) for 
C Oph among other OB stars in their sample. They found 
that in ( Oph the DEM sharply peaks at about 6 MK. While 
this is a significantly lower temperature than found for the 
DEM peak in case of 1 OriC (50 MK), it is higher than 
found for other stars of similar spectral types (^3MK). 
Thus, considering the X-ray temperature of COph, it is 
not straightforward to attribute its X-ray emission to the 
MCWS. On the other hand, recent studies of O stars with de- 
tected magnetic fields (e.g., HD 191612, HD 108) show that 
their X-ray properties are diverse (Naze et al. 120041 [2010) 
and may be difficult to fully reconcile with the predictions 
of the MCWS model. 

Clearly, new observations are needed to better under- 
stand the magnetic field of C Oph and its link with the X-ray 
emission from this star. 
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Fig. 5 Fe XVII (upper panel) and Ne X (lower panel) lines 
observed in the spectrum of COph (co-added MEG ±1). 
Vertical dashed lines indicate the rest-frame wavelength, 
corrected for the radial velocity. 
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